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FOREWORD 


This  report  OSURF  number  2467-2,  was  prepared  by  the  Electro- 
Science  Laboratory,  Department  of  Electrical  Engineering,  The  Ohio 
State  University  at  Columbus,  Ohio.  Research  was  conducted  under 
Contract  F44620-67-C-0095  of  the  Air  Force  Cambridge  Research 
Laboratories  at  Bedford,  Massachusetts.  Dr.  John  K.  Schindler  was 
the  Program  Monitor  for  this  research. 


ABSTRACT 


The  transient  fields  produced  by  vertical  and  horizontal  electric 
dipoles  with  impulsive  current  excitation  in  the  presence  of  a  half- 
space  are  considered.  For  special  source  and  observer  configurations, 
exact,  closed  form  expressions  for  the  fields  produced  by  impulsively 
excited  dipoles  above  a  nondispersive  half- space  are  given.  For  the 
dispersive  case,  estimates  of  the  transient  fields  from  asymptotic 
approximations  are  given.  A  knowledge  of  the  interactions  of  dipole 
antennas  and  a  half- space  is  useful  for  studies  of  the  background  return 
from  the  half- space  and  for  studies  of  the  scattering  by  small  targets 
on  or  near  the  half- space. 
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THE  TRANSIENT  FIELDS  OF  DIPOLE  ANTENNAS  IN  THE 
PRESENCE  OF  A  DIELECTRIC  HALF-SPACE 


I.  INTRODUCTION 

When  a  radar  target  is  located  on  the  ground  under  a  vegetation 
canopy,  its  scattering  characteristics  can  be  greatly  altered  from  those 
of  the  same  target  in  free  space.  The  ecology  of  the  ground  and  vege¬ 
tation  constitute  a  single  entity  but  it  is  more  convenient  to  consider 
separately  the  effects  of  the  ground  and  vegetation  cover.  A  Rayleigh 
region  leaf  model  of  the  vegetation  canopy  has  been  developed,  *  and 
calculations  of  the  effective  constitutive  parameters  in  the  30  to  300  MHz 
range  obtained.  In  this  frequency  range  the  effects  of  the  ground  can 
often  be  estimated  by  treating  the  ground  as  a  smooth  dielectric  surface. 
It  is  appropriate  therefore  to  consider  the  fields  produced  by  dipole 
antennas  in  the  presence  of  a  half- space.  Since  the  dipoles  may  be 
used  to  represent  either  sources  or  small  scatterers,  a  knowledge  of 
the  interactions  of  the  dipole  with  the  half- space  yields  estimates  of 
both  the  background  return  from  the  ground  and  the  effects  of  the  ground 
on  the  scattering  characteristics  of  a  target. 

A  general  transient  excitation  of  the  dipoles  is  considered  because; 
(1)  a  knowledge  of  the  ground  return  for  a  general  transient  excitation 
may  prove  useful  in  deciding  what  specific  interrogating  waveform  or 
frequency  spectrum  will  minimize  the  problems  of  background  return, 
and  (2)  a  general  extension  of  the  time  domain  concepts  developed  for 
free  space  targets*  to  the  problem  of  targets  located  in  a  natural 
environment  is  being  studied. 

The  first  tv/o  sections  of  this  report  consider,  respectively,  the 
transient  fields  of  vertical  and  horizontal  dipoles  above  a  nondispersive 
half- space.  For  the  vertical  dipole,  exact,  closed  form  solutions  for 
the  transient  fields  when  the  dipole  is  impulsively  excited  are  given 
when  the  observer  is  on  the  vertical  axis  containing  the  source  dipole 
or  when  both  source  dipole  and  observer  are  on  the  interface.  These 
results  are  carried  out  from  an  anisotropic  half- space  where  the  medium 
has  a  uniaxial  permittivity.  For  the  impulsively  excited  horizontal  dipole 
exact,  closed  form  solutions  when  the  observer  is  on  the  vertical  axis 
containing  the  dipole  are  given.  For  simplicity,  these  results  are  re¬ 
stricted  to  an  isotropic  half- space  and  observation  points  above  the  half¬ 
space.  These  simplifications  are  not  necessary,  however,  and  the  re¬ 
sults  can  be  extended  to  an  anisotropic  half- space  and  observation  points 
in  the  lower  medium. 
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If  the  half- space  lossy  or  stratified,  exact  solutions  for  the 
transient  fields  are  not  possible.  The  third  section  of  this  report  con¬ 
siders  long  and  short  time  estimates  of  the  transient  fields  obtained 
from  quasi-static  and  asymptotic  approximations.  Both  analytical  esti¬ 
mates  and  a  Fourier  synthesis  technique*  are  used  to  illustrate  some  of 
the  effects  of  a  finite  conductivity  on  the  transient  waveforms. 

The  final  section  of  the  report  discusses  practical  applications  of 
the  theoretical  results  in  predicting  the  effects  of  the  ground  on  the  scat¬ 
tering  characteristics  of  targets  located  on  the  ground.  A  theoretical 
model  for  the  scattering  by  targets  located  on  or  above  the  ground, 
including  interactions,  is  described.  A  linear  difference  equation  for 
the  transient  waveforms  is  illustrated,  and  the  application  of  this 
difference  equation  in  a  detection  and  discrimination  scheme  is  shown. 
Research  on  the  problem  is  continuing  on  Contract  FI 9628-67- C-0239i 
and  additional  results  will  be  reported  there. 


11.  VERTICAL  ELECTRIC  DIPOLE 

A  vertical  electric  dipole  located  at  a  height  h  above  an  anisotropic 
half-space  is  shown  in  Fig.  1.  It  is  assumed  that  the  half-space  z  >  0  is 
free  space  and  the  half- space  z  <  0  is  a  uniaxial  material  with  a  complex 
conductivity  tensor 

/o-y  0  0 

(1)  (ar)=:  0  (Th  0 

VO  0  (T 


where 

A 

^■v  =  +  j  «  «V 


The  subscripts  h  and  v  refer  to  the  horizontal  and  vertical  directions 
respectively.  The  problem  corresponds  to  that  considered  by  Wait* 
when  the  anisotropic  intermediate  layer  considered  by  Vvait  has  an 
infinite  depth.  It  can  be  shown  that*  for  the  special  case  of  a  nondis- 
persive  half-space,  i.e.,  a  half-space  which  can  be  characterized  by 
frequency  independent  indices  of  refraction  n^  and  ny,  and  an  observa¬ 
tion  point  located  on  the  z  axis  (directly  above  or  below  the  source 
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Fig.  !•  Dipole  over  anise  tropic  half- space. 


dipole)  that  the  real  time -dependent  electric  field  corresponding  to  an 
excitation  of  the  source  dipole  of  length  ds  with  a  current  It>(t)  is  given 
by 


E®  (z,  t)  =  — i _ Ilbll 

4TT<o(z+b)  c(z+h)  njj+1 


z  >  0 


where  the  dipole  is  of  length  ds  and  the  quantity  D(z,  t)  is  given  by 


D(Z|  t)  - 


_  2nj^n^ 


(z+h)2B*  [  A3 

•»  A/1  nn2\  2T 
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B  =  nynhT  +  A  ,  and 
F  -  nynhT  -  A 


In  Fq«  (2)»  the  direct  term*  ue.»  the  term  present  in  the  absence  of 
the  hal£> space,  has  been  removed  and  the  observation  point  is  re¬ 
stricted  to  the  half- space  i  >  0.  The  refractive  indices  nh  and  ny  are 
defined  respectively  by  the  ratio  of  horizontal  or  vertical  propagation 
constants  to  the  free  space  propagation  constant*  The  scattered  or 
interaction  field  in  £q.  (2)  is  exact  for  a  nondispersive  half- space  and 
an  observation  point  on  the  vertical  axis . 


The  initial  and  final  values  of  the  expression  in  £q.  (3)  are  easily 
obtained  as 

(4)  f  4nh(n^^-  I) 

V  c  I  (z+h)*  1  [nv(nh+‘)]‘  "h+l 


and 
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D(z,  ®)  = 
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D(z«t)  in  £q.  (3)  is  a  monotonically  decreasing  function  with  a  time 
constant  reoughly  equal  to  z+h/c.  Thus  while  the  expression  in  Eq.  (Z) 
is  not  simple,  the  resultant  waveform  is  a  relatively  simple  function 
with  an  initial  singularity  and  jump  discontinuity  and  then  a  smooth 
decay  to  the  final  vcilue*  A  sketch  of  the  waveform  of  Eq.  (2)  is  shown 
in  Fig.  2. 

For  the  special  case  of  a  perfectly  conducting  lower  half-space 

iny  =  njj  —  oo) 


which  agrees  with  an  image  theory  solution. 
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Fig.  2.  Sketch  of  transient  field  of  impulsively  excited 
vertical  electrical  dipole  above  half- space. 

Observation  point  on  vertical  axis  above  half¬ 
space. 

The  background  return  from  the  anisotropic  half- space  is  given 
by  setting  z  =  h  in  Eq.  (2),  Note  that  this  result  may  also  be  viewed  as 
the  interaction  of  a  vertical  electric  dipole  at  a  height,  h,  and  an  aniso¬ 
tropic  half- space. 

The  general  method  used  to  obtain  the  result  in  Eq.  (2)  can  also 
be  applied  for  observation  points  on  the  vertical  axis  in  the  lower  half¬ 
space  (z  <  0).^  In  this  case  an  isotropic  lower  half-space  (n^  =  nj^  =  n) 
is  considered  although  this  simplification  is  not  necessary  to  obtain  an 
exact,  closed  form  solution  for  the  time -dependent  field.  It  can  be 
shown  that  the  time -dependent  Hertz  vector  in  the  lower  medium  for 
observation  points  on  the  vertical  axis  is  given  by 
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(7) 


z*  +  h* 


+  zhj 


2  ^ni+1 

(ct)^ 


!!^;j 

•  [z  +  n*  h  +  (h+n*  z)J  1  +  (n*  -  1  )(h*  -  z*)/  (ct)*^  *  ‘  u  j  . 


The  tizne-^Jependent  electric  field  given  by 

(8)  E,z{z.t)  =  -  g]  *  ly  n,^{z.  t)  . 

is  quite  complicated,  as  is  obvious  from  Eq«  (7),  and  will  not  be  given 
here*  An  exact,  closed  form  expression  for  this  field,  however,  can 
be  obtained  from  the  results  in  Eqs.  (7)  and  (8)*  The  final  value  of  the 
field  is  given  by 


(9) 


Eia(2,co) 
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A  sketch  of  the  Hertz  vector  given  in  Eq.  (7)  is  shown  in  Fig.  3.  Note 
that  this  waveform  is  monotonic  decreasing  and  relatively  simple  de- 
I  pite  the  complexity  of  Eq.  (7).  It  follows  that  the  time -dependent  field 
given  by  Eq.  (8)  will  have  singularities  and  a  discontinuity  at  the  wave¬ 
form  origin  but  will  also  be  monotonic  decreasing  for  longer  times. 

One  other  special  orientation  of  source  dipole  and  observer  yields 
an  exact,  closed  form  expression  for  the  time -dependent  field.  If  both 
the  source  dipole  and  the  observation  point  lie  on  the  interface  (z  =  h  =  0), 
then  it  can  be  shown^  that 


where  p  is  the  horizontal  separation  of  the  source  dipole  and  observation 
point*  A  corresponding  expression  for  the  <|>  component  of  the  field 
(p,  <>,  z  cylindrical  coordinates)  is 


the  p  component  of  the  electric  field  cannot  be  obtained  from  this  pro 
cedure  because  the  z  dependence  of  the  intermediate  Hertz  vector  is 
unknown*  Van  der  Pol*  has  obtained  similar  results  for  the  case  of 
an  isotropic  half- space. 
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Fig.  3.  Sketch  of  t^me- dependent  Hertz  vector  of  impulsively 
excited  vertical  electric  dipole*  vertical  axis.  Source 
dipole  above  half- space*  observation  point  within 
half- space. 
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The  results  in  both  £q>  (10)  and  £q.  (11)  show  a  double  wavefront* 
i.e.,  one  wave  traveling  with  a  velocity  c  along  the  boundary  and  a 
second  wave  traveling  with  a  velocity  c/n^.  along  the  boundary.  7 he  lead¬ 
ing  terms  in  both  Eqs.  (10)  and  (11)  have  typical  p"^  behaviour  for  sur¬ 
face  waves.  Note  that  the  ^  component  of  the  field  (Eq.  (H))  has  a  finite 
duration.  A  sketch  of  the  waveform  in  Eq.  (11)  is  shown  in  Fig.  4.  The 
time -dependent  fields  given  by  Eqs.  (10)  and  (11)  are  typical  of  the  excita¬ 
tion  produced  by  a  vertically  polarized  source  at  very  low  grazing  angles. 
The  results  in  Eqs.  (10)  and  (11)  are  for  the  total  fields,  i.e.,  the  direct 
term  is  included. 


III.  HORIZONTAL  ELECTRIC  DIPOLE 

if  the  vertical  electric  dipole  in  Fig.  1  is  replaced  with  a  hori¬ 
zontal  electric  dipole,  oriented  in  the  x  direction,  the  problem  has  an 
added  complexity  because  of  the  loss  of  a  <t>  symmetry.  It  is  convenient. 


8 


Fig-  4.  Sketch  o£  transient  field  of  impulsively  excited 
vertical  electric  dipole*  Source  dipole  and 
observation  point  on  interface- 


but  unnecessary,  to  assmne  an  isotropic  lower  half- space*  For  a  dis¬ 
persionless  half- space  and  an  observation  point  on  the  vertical  axis  in 
the  upper  medium  (z  ^  0)>  it  can  be  shown  that  for  a  dipole  of  length  ds 
the  time -dependent  field  produced  by  a  current  excitation  l^t)  is  given 
by* 


(12) 


z  >  0 
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c(z+h)  ^  n+1  (z+h)* 


The  final  value  of  the  magnetic  field  is  zero  and  the  final  value  of  the 
electric  field  is  given  by 


(14) 


E^^(z.<») 


Ids  n^  -  1 

4Tr€o(z+h)^  t?TT 
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A  sketch  of  the  waveform  given  in  Eq.  (13)  is  shown  in  Fig#  5.  For  the 
results  in  Eqs.  (12)  and  (13),  the  direct  term  has  been  removed.  Similar 
results  under  the  same  restrictions  for  the  time-dependent  fields  in  the 
half- space  z  <  0  can  be  obtained,  and  the  results  for  the  horizontal  dipole 
can  also  be  extended  to  an  anisotropic  half- space. 

Details  of  the  methods  used  to  obtain  the  results  for  the  vertical  and 
horizontal  dipoles  are  given  in  Reference  4  and  will  be  included  in  a 
future  report.  Basically,  the  inverse  Fourier  transform  of  the  Sommer- 
feld-type  integral  solutions  is  taken  after  a  change  of  variable.  It  turns 
out  that  for  a  nondispersive  half-space,  the  inverse  transform  may  be 
obtained  for  vertical  or  horizontal,  electric  or  magnetic  dipoles,  for 
observation  points  in  either  half-space,  assuming  an  interchange  of  orders 
of  Integration  is  valid.  For  arbitrary  locations  of  the  source  dipole  and 
observer,  the  method  results  in  a  remaining  finite  integral  which  must 

be  evaluated  numerically.  For  the  special  cases  considered  here,  the 
remaining  integral  can  be  evaluated  exactly  in  closed  form. 
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Fig,  5,  Sketch  of  transient  field  of  impulsively  excited 
horizontal  electric  dipole  above  half- space. 
Observation  point  on  vertical  axis  above 
half- space# 


The  most  significant  difference  between  the  interaction  of  vertical 
and  horizontal  electric  dipoles  and  a  nondispersive  half- space  is  the 
addition  of  a  doublet  termi  which  varies  as  the  inverse  of  twice  the 
dipole  height  above  the  half-space*  for  the  horizontal  dipole. 


SIMPLIFIED  ESTIMATE  OF  RESULT 
FOR  VERTICAL  DTPOLE 


It  is  evident  from  the  complexity  of  the  results  for  the  time-depend¬ 
ent  fields  and  Hertz  vectors  presented  in  the  first  two  sections  that  an 
analytical  utilization  of  the  expressions  would  be  extremely  cumbersome. 
However,  sketches  of  the  waveforms  deduced  from  an  examination  of  the 
expressions  are  quite  simple.  This  suggests  that  with  the  singularities 
in  the  waveforms  removed,  very  simple  approximate  expressions  for  the 
waveform  may  be  feasible.  For  a  homogeneous  half- space,  the  transient 
fields  for  either  vertical  or  horizontal  dipoles  must  be  quite  simple  in 
form  once  the  singularities  associated  with  the  interface  are  removed. 


Using  the  known  initial  and  final  values  of  the  waveform  for  the 
vertical  electric  dipole  in  Eq.  (3),  an  approximation  of  the  form 

(15)  D(z,t)a<A+Be 


where 

A  =  _L_  r^»h<V-^)  ^  1 

(a+h)'‘  [(nv(nh+l))*  J 


is  suggested.  In  Eq.  (15)  the  delay  must  be  (z+h)/c,  the  time  constant 
could  in  general  be  different  but  is  set  to  (z+h)/ c  here  from  observations 
of  the  exact  waveform. 

To  test  the  approximation,  the  waveforms  predicted  by  Eqs.  (3)  and  (15) 
arc  compared  in  Figs.  6  and  7  for  various  refractive  indices  and  both 
isotropic  and  anisotropic  half-spaces.  It  is  clear  from  these  results 
that  Eq.,  (15)  is  a  reasonable  approximation  for  refractive  indices  some¬ 
what  greater  than  ^2  for  both  isotropic  and  anisotropic  half-spaces.  A 
transformation  of  Eq.  (2)  using  Eq.  (15)  to  replace  Eq.  (3)  leads  to  a 
simple  expression  for  the  frequency-dependent  field  below  a  vertical 
dipole  which  should  be  a  reasonable  approximation  over  the  entire 
spectrum.  Admittedly,  in  most  cases  one  is  rarely  inter^-.sted  in  the 
fields  directly  below  a  vertical  dipole;  the  value  of  these  results  in  that 
sense  lies  in  the  demonstrated  feasibility  of  using  a  simple  time  function 
to  replace  the  complicated  exact  expression.  They  suggest  that  from  a 
knowledge  of  the  singularities  and  the  initial  and  final  values  of  the  time- 
dependent  waveform,  simple  estimates  of  the  waveforms  for  more  inter¬ 
esting  observer  locations  may  be  possible.  It  will  be  shown  later  that 
the  dipole -ground  interactions  obtained  in  the  first  two  sections  are  two 
of  the  four  interaction  results  needed  for  a  theoretical  model  of  the  scat¬ 
tering  from  small  ta  '  ?  located  on  the  ground.  Thus  for  the  problem 
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(2>hr  Oth.z.t) 


g,  7.  Comparison  of  exact  and  approximate  transient  fields  of 
impulsively  excited  vertical  electric  dipole  above  half¬ 
space.  (a)  nj^  =  4,nv  =  2;  (b)  nj^  =  ny  =  3.46;  (c)  nj^  =  3, 
n^  “  2j  (d)  nj»  —  5|  n^  —  2 » 


at  hand,  the  transient  fields  on  a  vertical  axis  containing  the  dipole 
source,  and  approximations  for  them,  are  of  more  than  academic 
interest. 


V.  EFFECTS  OF  LOSS 

If  the  half- space  is  dispersive  or  stratified,  the  techniques 
utilized  in  the  previous  sections  are  not  useful  for  obtaining  the  transient 
fields.  In  these  cases,  even  when  the  dielectric  constant  and  conductivity 
of  the  half- space  are  assumed  to  be  constant  with  frequency,  an  exact 
closed  form  solution  for  the  transient  fields  does  not  appear  to  be  possible. 
An  exception  is  noted  for  a  particular  dispersive  half- space,  ^  rather 
unrelated  to  the  ground  effects  of  interest  ► 

If  the  half- space  dielectric  constant  and  conductivity  are  assumed 
constant  with  frequency  and  attention  is  confined  to  sources  and  observa¬ 
tion  points  above  or  on  the  half- space,  then  :ertain  deductions  concern¬ 
ing  the  form  of  the  transient  fields  are  possible.  In  the  limit  of  long 
times  (low  frequencies)  the  fields  must  decay  to  their  static  values,  and 
the  long  time  values  of  the  transient  waveforms  can  be  obtained  from 
quasi- static  estimates  of  the  frequency -dependent  fields  and  the  final 
value  theorem.  A  more  sophisticated  long  time  estimate  (holding  for 
shorter  times)  is  obtained  by  assuming  the  displacement  currents  axe 
negligible.  That  is,  the  propagation  constant  of  the  half- space  is 
approximated  as  kj  =  -  ,  This  approximation  has  been  used  in 

most  of  the  analytical  solutions  for  the  transient  fields  of  dipole  sources 
in  infinite  and  semi-infinite  conducting  media. ^  For  short  times, 
estimates  of  the  fields  from  high  frequency  asymptotic  approximations 
are  indicated.  If  the  observation  point  and  source  dipole  are  confined 
to  be  above  or  on  the  half- space  and  only  the  secondary  (fields  due  to 
the  presence  of  the  half- space)  transient  fields  are  considered,  then 
these  transient  fields  must  be  characterized  by  singularities  and  a  jump 
discontinuity  at  the  initial  time  of  arrival  of  the  field  followed  by  a  smooth 
decay  to  the  final  static  value.  That  is,  in  general  when  the  waveform 
first  departs  from  zero  singularities  of  various  order  and  a  jump  dis¬ 
continuity  will  occur,  but  beyond  this  initial  arrival  time  the  geometry 
of  the  problem  is  such  that  only  a  smooth  decay  of  the  waveform  is  poss¬ 
ible.  The  above  assumes  the  surface  type  fields,  i.e.,  fields  arriving 
partially  via  travel  within  the  half- space  can  be  neglected.  For  the 
asymptotic  approximations  of  interest  this  means  roughly  that 
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where  Cr  is  the  relative  dielectric  constant  of  the  half* space. 

In  this  report,  attention  is  confined  to  the  final  values  of  the 
transient  fields  and  to  the  short  time  estimates  predicted  by  high  fre¬ 
quency  asymptotic  approximations.  The  asymptotic  approximations 
utilized  are  those  given  by  Banos.* 


Lone- Time  Estimates 


Bannister*  has  given  quasi -static  estimates  of  the  frequency- 
dependent  fields  of  vertical  and  horizontal  dipoles  above  a  conducting 
half-space.  The  geometry  of  the  general  problem  is  shown  in  Fig.  8. 
The  long  time  values  of  the  total  transient  fields  for  impulsive  current 
excitation  for  vertical  and  horizontal  electric  dipoles  at  an  arbitrary 
location  above  or  on  the  half-space  are  given  below. 


Fig.  8.  Geometry  of  general  problem. 
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Vertical  Electric  Dipole  (h  >  0,  a  >  0) 


Ez(p.  z.“)  =  - 


|s_  \  1  A  3(z-h)^'\  1  A  3(z-fh)*  VI 


where 


Ro*  =  p*  +  (z  h)' 


Ri*  =  p*  +  (z+h)^ 


ri|^±)  +  3£U 

•  ^TTCo  L  R: 


J 


Hj.(p,  z,  c»)  =  0 


Horizontal  Electric  Dipole  (p  >  z+h,  z  >  0,  h  >  0) 


piP.<P.  )  4^,^  1^^  ^ 


E4,(p.  <)> 


»  IdF  sin(|> 
,  z,oo)  =  — ^ - 1. 

4ireo 


[s? 


Hp(p,c;>,  Z,oo)  =  H  {p,(|),  Z,cc)  =  H2(p,9»  Z,oo)  =  0 


Note  that  the  above  expressions  contain  both  the  primary  and  scattered 
fields*  A  comparison  of  the  long  time  values  obtained  here  with  those 
for  a  nondispersive  half-space  illustrates  that  if  the  half-space  has 


some  conductivity  the  final  values  of  the  transients  are  independent  of 
the  half- space  parameters*  Similar  results,  i*e«,  long  time  solutions, 
for  vertical  and  horizontal  magnetic  dipoles  can  be  obtained* 


B.  Short  Time  Estimates 

Asymptotic  approximations  of  the  integral  expressions  resulting 
from  an  exact  solution  of  the  dipole,  half- space  problem  are  somewhat 
more  difficult  to  obtain  than  those  for  the  quasi -static  limit*  One  also 
finds  that  such  approximations  are  usually  specialized  to  various  regions* 
Banos,  *  for  example,  obtained  three  term  asymptotic  approximations  in 
the  vicinity  of  the  vertical  axis  and  in  the  vicinity  of  the  surface  inter¬ 
face  but  only  two  term  expansions  over  the  entire  upper  hemisphere* 


1 .  Vertical  axis 

Banos  gives  asymptotic  approximations  in  the  vicinity  of  the 
vertical  axis  which  reduce  to  the  following  forms  for  observation  points 
on  the  vertical  axis  (p  =  0)  and  on  or  above  (z  >  0)  the  conducting  half¬ 
space* 


Vertical  Electric  Dipole  (p  =  0,  z  ^  0) 


(23) 


e|{0,  z) 


Zplds 

Ztt 


-j-(h+z) 

c 


6 

-1  +  3n  -  5n’  +  3n^ 

|+Ji2Ei_  . 

kj  -(h+z)* 

1+n 

'  (1)  k,{h+z)*  (1+n)^ 

where  is  the  free  space  impedance, 

k. 


n  - 

n  -  —  w 


*  Note  change  in  definition  from  that  used  in  first  two  sections. 


El  =  -  2  +  4n  -  5n*  +  3n*  . 


and  €y.  O'  are  the  relative  dielectric  constant  and  conductivity  of  the 
half- space  respectively-  The  above  result  is  for 


!(?)■  «r  -  j“|Ao<fi(z+l')  I  >1-0 


For  a  dielectric  half- space  (oi=  0)  Eq.  (23)  yields  a  transient  waveform 
(25) 


E^(a.t')  =  — 

*  2Tr(h+z)* 


tllL  6(t')  + 


6c*t’ 


4^+1 

/I.  5 


c  f  3-4/  4^^ 

i^V  4^+1  / 


^3-4^1 


6c\'^ 

j€7(h+z)^ 


/  3  5  \  *1 

u(t')  . 

where  t'  =  t  -  (z+h/c)- 

It  can  be  shown  that  the  latter  result  agrees  with  Eq.  (2)  for  t‘  equal 
zero.  That  is,  for  a  dielectric  half- space  both  the  delta  function 
singularity  and  jump  discontinuity  are  correctly  predicted  by  the 
asymptotic  approximation. 


Horizontal  Electric  Dipole 

For  the  horizontal  electric  dipole,  Banos  gives  an  asymptotic 
expansion  for  observation  points  above  or  on  the  half-space  and  in  the 
vicinity  of  the  v^ertical  axis  for  the  condition  in  Eq.  (24),  On  the  vertical 
axis,  this  result  reduces  to 
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E^U)  = 


Zq  Ids  I  c 


2n  |2(z+h) 


fe)*  (^) 


1 _  r2n*  -  6n^-f  2n*-h  n+  l1  c 

2w  j  I  1  l-n  J 

j  —iz+hr  ^ 


For  a  dielectric  half-space  (cri  =  0),  Eq.  (26)  yields  a  transient  waveform 


'Vo  Zq  Ids 

“  2ir(z+h) 


l/Ell] 

2c  W+l  .j 


2(z+h)*(J^+l) 


«r*  +  «r*^*  +  2ir  -  6JI7  +  2 


Note  that  the  first  two  terms  in  Eq.  (27),  i.e.,  the  doublet  and  impulse 
weights,  agree  with  those  in  Eq.  (12).  However,  the  coefficient  of  the 
step  discontinuity  in  Eq.  (27)  does  not  agree  with  that  obtained  from 
Eq.  (12)^^  v/hen  the  differentiations  are  carried  out.  It  is  shown  in 
Reference  11  that  the  coefficient  from  Eq.  (12)  is 


- ; - r—= - +  ,  _ - [.itr*+ 6er*/*+5*r  +  J^  +  2r 

2{z+h)*2ff  ^  4  ^ 


If  the  half- space  is  dispersive,  the  long-time  values  of  the  transient 
fields  given  in  Eqs.  (16)  -  (22)  hold  for  any  observation  point  in  the  upper 
hemisphere  z  >  0.  For  observation  points  on  the  vertical  aixis  in  the 
upper  hemisphere,  some  properties  of  the  short-time  behaviour  can  be 
obtained  from  the  asymptotic  results  in  Eqs.  (23)  and  (26).  For  the 
vertical  electric  dipole,  the  weight  of  the  impulse  singularity  at  t'  equal 
to  zero  is 


f,(z,  0)  = 

2iT(h+z) 


Ed. 


and  the  step  discontinuity  is 


(30) 


fj(z.O) 


ZqMs 

2Tr{h+z)*(l+JTr') 


m 

i  • 


(T _ 


h+2 


+  ^ 


For  the  hc^  ii&ontal  electric  dipole,  similar  utilization  of  the  initial  value 
theorem  yields 


(31) 


fo(2.  0) 


S^Qlds  1 

- -  -  • 

2c(z+h) 


> 


for  the  weight  of  the  doublet  term  and 


(32) 


f,(z,  0)  = 

2ir 


_,^+l 


for  the  weight  of  the  impulse  term.  Note  that  for  both  the  vertical  and 
horizontal  d:poles,  the  weight  of  the  highest  order  singularity  is  identical 
to  that  for  the  corresponding  dielectric  half- space  regardless  of  the 
conductivity  of  the  half-space.  For  the  horizontal  dipole,  the  magnitude 
of  the  jump  discontinuity  is  not  obtained. 

To  establish  an  estimate  of  the  transient  waveforms  predicted 
by  the  asymptotic  estimates  in  Eqs  •  {Z3)  and  (26)  a  Fourier  synthesis 
technique,  ^  i.e.,  a  numerical  inversion  of  the  frequency-dependent 
expressions  using  calculations  at  harmonically  related  frequencies,  can 
be  used.  In  using  this  technique,  the  period  of  the  fundcunental  frequency 
IS  selected  to  be  at  least  10  times  the  relaxation  time  (Tr  =  c/cr),  and  at 
the  same  time  the  relative  paraiueters  of  the  half-  space  and  the  dipole 
height  are  selected  to  satisfy  Eq.  (24).  In  Fig.  9,  the  normalized  back¬ 
ground  or  backscattered  response  waveforms  for  a  vertical  electric 
dipole  for  a  height  z  =  h  =  200  meters,  half- space  relative  dielectric 
constant  of  12  and  conductivities  of  0.1,  0.01  and  0,001  mhos  per  meter 
are  shown.  The  abscissa  scale  is  in  units  of  the  relaxation  time  of  the 
half- space,  and  the  coefficient  of  the  impulse  term  (Eq.  (29))  has  been 
subtracted  from  the  calculations  before  synthesis.  Corresponding  wave¬ 
forms  for  the  horizontal  dipole  for  the  same  parameters  are  shown  in 
Fig.  10,  Note  that  in  this  case  the  coefficient  of  the  doublet  term  (Eq.  (31)) 
has  been  subtracted  and  also  a  reflection  coefficient  type  expression 
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2-/W  h  e;  {h,  »/t,) 


-I  0  I  2  ?  4  5 

UNITS  OF  RELAXATION  TIME  Tr  ^ 


Fig.  9.  Transient  backscattered  fields  of  impulsively  excited 
vertical  electric  dipole  "above  conducting  half- space. 
€j.=  12,  h=200  m,  cr  =  0.1,  0.01,  0.001.  Impulse 
singularity  remDved. 
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h  =  200M 

€r=l2.0 

-  o-=  0.001 

-  cr  =  0.0l 

••  <T=0.l 


UNITS  OF  RELAXATION  TIME  — 

Fig.  10,  Transient  backscattered  fields  of  impulsively  excited 
hori/.ontal  electric  dipole  above  conducting  half- space. 
€j,=12,  h  =  200m,  (r=0.1,  0.01,  0.001.  Doublet 
singularity  and  portion  of  impulse  singularity  removed 


(33) 


has  been  subtracted  before  synthesis.  It  is  evident  from  Fig.  10  that 
the  impulse  singularity  weight  is  not  independent  of  the  half- space 
conductivity  since  the  waveforms  are  still  dominated  by  an  impulse  type 
singularity.  An  estimate  of  the  magnitude  of  vbe  jump  discontinuity 
could  be  obtained  by  subtracting  both  Eqs.  (31]  and  (32)  before  synthesis. 
However,  it  is  clear  from  Fig.  9  and  Eq.  (30)  t  'at  the  synthesis  pro¬ 
cedure  cannot  faithfully  estimate  the  discontinuities  involved  for  con¬ 
ductivities  in  this  range  except  possibly  by  using  many  more  harmonics 
than  the  40  used  for  the  results  in  Figs.  9  and  10,  For  either  the  vertical 
or  horizontal  dipoles,  the  essential  form  of  the  transient  waveforms  is 
clear  but  an  extensive  computational  program  would  be  required  to 
establish  exact  details  of  the  waveforms  over  a  range  of  half-space 
parameters,  source  dipole  heights  and  observer  locations.  Summarized 
below  are  the  established  initial  and  final  values  of  the  transient  wave¬ 
forms  for  both  vertical  and  horizontal  dipoles  for  observation  points  on 
the  vertical  axes  and  in  the  upper  hemisphere  (z  >  0),  These  values  are 
for  the  z  component  of  the  electric  field  for  the  vertical  electric  dipole 
and  horizontal  component  of  the  electric  field  parallel  to  the  dipole 
orientation  for  the  horizontal  dipole. 


Vertical  Electric  Dipole 


Half-  space 


Initial  Value  (jump  discontinuity)  Final  Value 


Perfect  Conductor 


Lossless  Dielectric 


Conducting 


Ids 

Ids 

2TT€o(z+h)^ 

2TT€Q(z+h)^ 

Ids 

Ids 

27reo(z+h)* 

2  TT€q(z +}!)■* 

Cy+i 

Ids 

2TT€Q(h+z)^ 

fr(.Fr+l)' 

r 

+ 


h+z 


-I- 


r 


_ Ids 

2ir€oT^+hp 
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Horizontal  Electric  Dipole 


Initial  Value 


Final  Value 


Perfect  Conductor 

4ir€o(z+h)" 


Ids 

4Tlto(2+h)^ 


Lossless  Dielectric 


Ids  c 


Conducting 


2Tr€o2(z+hr|j€j(4r;+l) 
4rr^(J??+i)*  * 

♦  [2«J.*  +  6<r*/  *  +  5*r  +  ^ + 2] 


Ids  f^r  •  ^ 
4Tr€Q(z+h)^  I  €j.+l 


4Tr€o(z-!-h)^ 


2  •  Entire  upper  hemisphere  (z  >  C) 

•v  - 

Banos  gives  two  term  asymptotic  expansions  for  observation 
points  anywhere  in  the  upper  hemisphere  z  >  0, 


Vertical  Electric  Dipole 


^ki“'  -  —  't;  *  rT'  - '  -  w — 


3(;h+z)*  3(h+z)^ 

- ~  ^ 


+  2  cos  0 _  .  ^o  ^  )  3(h-i  z)^ 

cos  6  4  n>!l-n^sin^0  ^  Rl  Ri^  Kq  Ri^  R^ 


.  3(h4-z)^  E,  J  >j  {koVl'?')^+2)  5(h4z)^ 

‘  ^  Rr  ■  VRi^  ^ 


(34) 

(cont>) 


where 


8(h+z)* 


r  ■■; ,  jcoto  — 


u(e  -  0q) 


>  » 


h+z  =  P.icos0  ,  p  =  Risin0 


®o  =f-  +  |n|  +1  |n|*+  o|nh 


El  =  -  1  +  !I 


sin^O  |n*cos0+  4^  -  sin^O  1  ^ 


l-n^sin^0  Lv:os  0  +  n  Jl  -  sin^ 0  J 


cos  0  +  ^  f  3  sm^0  A  ^  sin^0  cos^0  sin^ 0 

^  ^  \  ^  (1  -  sin^  B)^^ 

cos  0  +  n  ^  1  -  s  In^  0 


2  cos  0 


For  a  dielectric  half-space,  the  asymptotic  approximation  in  Eq,  (34) 
yields  a  transient  waveform 


(35) 


E|(e.  t) 


Ids  Zo 
4Tr 


Hh+zV  ^  2  trcose  / 1 h+zV  N  ‘ 

Ri  /  «j.cos0+  «r-sin^0\\^*  /  J  . 


Zeycos  0 

€3.  cos  0  +  Jcj.  -  sin^0 
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(35) 
(cont. ) 


+ 


2  u(e  -  Op) _ 


where  in  Ei  and  0q,  n  -**1/  The  last  term  in  Eq.  (34)  is  inverse 
transformed  under  the  condition  that  the  time  function  be  real.  This  pro¬ 
cedure  properly  yields  the  delta  function  singularity  at  t  =4€7p/c  when 
h  =  z  =  0  as  is  evident  from  Eq.  (34)  if  we  set  h  =  z  =  0  and  B  =  tt/2  before 
transforming.  The  waveform  in  Eq.  (35)  specialized  either  to  the 
vertical  axis  (p  =  0)  or  to  the  case  of  both  source  and  observation  point 
on  the  interface  (h  =  z  =  0)  correctly  predicts  the  singularities  known 
from  the  exact  solutions,  but  yields  incorrect  magnitudes  for  the  jump 
discontinuities. 

It  is  interesting  to  observe  the  precursor  for  the  surface  field 
(last  term  in  Eq.  (35))  when  either  the  dipole  source  or  observation 
point  is  raised  above  the  interface  and  0  >  0^.  Similar  effects  have  been 
noted^®  from  asymptotic  scattering  approximations  when  the  assumed  path 
is  not  a  minimum  time  path.  In  this  case,  the  delay  time  c  is  the 

maximum  delay  possible  for  a  v/ave  traveling  via  the  dielectric.  When 
shorter  delays  are  possible,  i.e.,  z  or  h  0,  the  asymptotic  approxi¬ 
mation  properly  predicts  a  surface-type  field  prior  to  JI^p/c.  Previous 
expedience  indicates  that  such  an  estimate  is  good  only  in  the  vicinity  of 
this  time  however. 


27 


Horizontal  Electric  Dipole 


^  For  the  horizontal  electric  dipole,  the  two  term  expansion  by 
Banos*  yields  * 


E®  (W  =  Jr  ]l£  +  *  +  3p2 

4ir  R,  ^  jkoR:*  r7 

jZnkoCOsO  ^ 
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Ip"  . 
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For  a  dielectric  half-spacet  the  asymptotic  estimate  in  Eq«  (36)  yields 
a  transient  waveform 


(37) 
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(37) 

(cont.) 


•  (t  -  T  )u(t  -  t)  dT 


J 

where  in  Ei',  Ei  and  0^,  n  —  l/^/€^  . 

Again,  specializing  Eq*  (37)  to  the  vertical  axis  (0=0,  p  =0)  one  finds 
that  the  singularities  agree  with  those  in  Eq.  (27)  but  the  jump  magnitude 
does  not  agree  with  the  result  given  in  the  Table.  We  conclude  that  the 
two  term  asymptotic  expansions  given  by  Banos  correctly  predict  only 
the  highest  order  singularities  in  the  transient  v/aveforms  for  a  dielectric 
half- space.  It  follows  that  the  dependence  of  the  jump  discontinuity  on 
the  conductivity  of  the  half- space  for  arbitrary  observation  points  cannot 
be  established  from  these  results.  A  numerical  inversion  of  Eqs ,  (34) 
and  (36)  with  the  singularities  removed  would  be  pointless  since  the 
results  for  a  dielectric  half- space  are  incorrect.  Formally,  one  could 
derive  a  third  term  in  the  expansion  using  the  procedure  given  by  Banos. 
Assuming  this  led  to  the  correct  jump  discontinuities  for  a  dielectric 
half- space,  dependence  of  the  jump  discontinuity  on  the  conductivity  might 
then  be  inferred  with  some  confidence.  For  observation  points  where  the 
last  three  terms  of  Eq.  (36)  are  needed,  i.e.,  0  >  Gq,  the  precursor  seen 
in  the  case  of  the  vertical  dipole  is  again  evident  except  when  both  the 
dipole  and  observation  point  lie  on  the  interface. 

For  both  the  vertical  and  horizontal  electric  dipoles,  application 
of  the  initial  value  theorem  to  Eqs.  (34)  or  (36)  yields  the  same  weights 
for  the  doublet  and  impulse  singularities  as  given  in  Eqs.  (35)  and  (37) 
respectively.  That  is,  the  two  term  asymptotic  expansions  predict 
doublet  and  impulse  singularities  independent  of  the  half- space  con¬ 
ductivity.  We  have  already  seen  that  for  the  horizontal  dipole,  an  im¬ 
pulse  weight  dependent  on  the  half- space  conductivity  is  obtained  from 
a  three  term  expansion  on  the  vertical  axis.  As  noted,  on  the  vertical 
axis  the  two  term  expansions  correctly  predict  the  singularities  for  a 
dielectric  half-space.  However,  for  a  dispersive  half- space  the  two 
term  expansions  are  only  reliable  in  predicting  the  weight  of  the  doublet 
singularity  in  the  transient  waveform.  Unfortunately,  the  two  term 
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expansions  giv'en  by  Banos  appear  to  be  the  most  advanced  asymptotic 
development  available  for  observation  points  in  the  entire  upper  hemi¬ 
sphere.  It  seems  clear  that  additional  terms  are  needed  in  order  to 
estimate  the  effects  of  the  half-space  conductivity  on  the  lower  order 
singularities  and  jump  discontinuities  of  the  transient  waveform.  Note 
that  this  specifically  requires  that  the  displacement  currents  not  be 
neglected  because  it  is  the  short  time  character  of  the  waveform  which 
is  of  interest.  Thus  the  numerous  developments  where  displacement 
currents  are  neglected^  are  not  useful  in  resolving  the  singularities 
and  discontinuities. 


3.  Vicinity  of  interface 

Banos  has  also  given  three  term  asymptotic  expansions  for  the 
fields  produced  by  both  vertical  and  horizontal  electric  dipoles  located 
above  a  conducting  half- space  for  observation  points  in  the  vicinity  of 
the  interface.  That  is,  it  is  assumed  for  these  results  that  both  the 
observer  and  source  dipole  heights  (z  and  h  respectively)  are  much 
less  than  the  horizontal  range,  p,  of  the  observation  point.  A  major 
portion  of  the  book  is  devoted  to  this  development  and  the  vertical  axis 
results  utilized  earlier  were  intended  as  a  supplement  for  small  hori¬ 
zontal  ranges  where  the  interface  results  are  not  applicable. 

Explicit  expressions  for  the  fields  in  the  vicinity  of  the  interface, 
for  either  vertical  or  horizontal  dipole  sources,  are  quite  lengthy  and 
will  not  be  given  here.  Note  that  Banos  does  not  present  explicit  ex¬ 
pressions  for  the  fields,  his  results  are  asymptotic  expansions  of 
certain  basic  integrals  from  which  the  fields  are  obtained  using  standard 
formulas®  involving  primarily  second  derivaf’  with  respect  to  the 
observer  height  or  horizontal  range.  No  simplification  of  the  expressions 
occurs  if  either  the  source  dipole  or  observer  is  located  on  the  inter¬ 
face,  and  one  notes  the  same  precursor  effects  found  with  the  two  term 
upper  hemisphere  expansions.  However,  if  both  observer  and  source 
are  on  the  interface,  then  the  expressions  for  the  fields  are  greatly 
simplified.  Also,  for  this  case  we  have  an  exact  result  for  the  vertical 
dipole  which  offers  a  check  on  the  asymptotic  development  in  the  limit 
of  zero  conductivity. 

For  the  vertical  electric  dipole,  if  both  source  and  observer  are 
on  the  interface,  then 
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where  in  the  three  term  expansion  for  the  integral  related  to  the  phase 
delay  k^p/n,  certain  terms  have  been  deleted.  It  is  possible  to  do  this 
judiciously  in  the  case  because  the  exact  solution  for  the  transient  fields 
for  a  dielectric  half- space  is  knov/n.  For  a  dielectric  half- space  the 
result  in  Eq.  (38)  yields  a  waveform 


where  the  terms  with  a  1/p*  dependence  have  not  been  written  out#  It 
can  be  shown,  by  adding  the  direct  term  to  Eq#  (39),  that  the  singularity 
eind  jump  at  t  =  p/c  and  the  singularity  at  t  =  Ji7 p/c  agree  with  the  exact 
solution  in  Eq#  (10)  for  an  isotropic  half-space.  It  is  interesting  to  note 
that  no  deletions  were  necessary  to  obtain  agreement  at  t  =  p/c  and  only 
one  deletion  was  necessary  to  obtain  agreement  for  the  singularity  at 
t  =  J^p/c*  However,  it  was  not  possible  by  any  inclusion  or  omission  of 
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terms  to  obtain  agreement  for  the  jump  discontinuity  at  t  =  JcJTp/c.  In 
the  asymptotic  expansions  of  the  integrals  related  to  the  phase  delays 
k^p  and  k^p/n,  terms  to  order  l/p^  were  retained  in  each  case. 


For  the  horizontal  electric  dipole  with  both  source  and  observer 
on  the  interface  and  retaining  terms  only  to  1/p^ 


For  a  dielectric  half-spacci  Eq.  (40)  yields  a  waveform 


For  a  horizontal  electric  dipole,  if  both  source  and  observer  are 
on  the  interface  (h  =  z  -  0)  and  the  half- space  is  nondispersive  an  exact 
closed  form  solution  for  the  radial  component  of  the  transient  field  is 
easily^  obtained  as 
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where 


T  =  [(.,«,  (f)‘  ■  .,] 
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For  the  same  configuration  and  a  conducting  half -space ,  Wait^  has 
given  an  approximate  solution  for  the  component  of  the  radial  field  in  the 
xz  plane  =  0),  neglecting  displacement  currents,  as 
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The  results  in  Eqs .  (42)  and  (43)  represent  the  extremes,  respectively, 
of  neglecting  the  conduction  and  displacement  currents.  From  Eq.  (40), 
the  weight  of  the  delta  function  singularity  for  the  seconary  field  is 
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and  the  magnitude  of  the  jump  discontinuity  is 


(45) 


\(P» t) 

2 


Ids  Z^ 
4Tr 


0* 


2c 


«o(«r“^)p^ 


From  the  results  in  Eqs.  (43),  (44)  and  (45),  a  reasonable  estimate  of 
the  radial  field  for  an  impulsive  current  excitation  for  all  times  should 
be  feasible-  Waiti^  has  plotted  Eq.  (43)  (second  term)  and  this  result 
could  be  corrected  for  short  times  using  Eqs.  (44)  and  (45). 
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VI.  SCATTERING  MODEL  FOR  TARGETS 
LOCATED  ON  A  HALF-SJPACE 

Effect  of  Ground  on  tlie  Dipole 

Modes  of  Scatterers 

1  •  Dipole  inodes  of  scatterers  in  free  space  . 

At  low  frequencies,  i.e.,  when  a  target  is  electrically  small,  it 
is  possible  to  represent  a  scatterer  by  short  electric  and  magnetic 
dipoles.  A  maximum  of  three  mutually  orthogonal  electric  and  three 
mutually  orthogonal  magnetic  dipoles  is  required  for  any  object,  and 
these  numbers  are  reduced  for  some  shapes.  The  magnitude  of  any 
one  dipole  moment  is  given  by  the  product  of  the  parallel  component  of 
the  incident  field,  Eq  (for  electric)  or  Hq  (for  magnetic),  and  the  polari¬ 
zability  of  the  target,  a.  At  very  low  frequencies,  the  polarizabilities  are 
constants  and  the  resultant  scattered  field  has  an  s^(s  =  p+jw)  Rayleigh 
dependence.  However,  if  the  polarizabilities  are  made  frequency  de¬ 
pendent  the  resultant  scattered  field  ic  valid  into  the  first  resonance  region 
as  well  as  the  Rayleigh  region.  For  example,  the  conducting  sphere  of 
radius,  a,  has  an  electric  polarizability,  and  a  magnetic  polarizability, 
CLj^»  given  by 


V(s)=-^ 

s+£ 

a 

where  and  are  constants,  and  c  is  the  speed  of  light. 

Consequently,  the  response  of  the  sphere  has  s-plane  poles  at 
(“  1/2  +  jN/3/2)c/a  due  to  the  electric  dipole  and  at  -c/a  due  to  the  mag¬ 
netic  dipole.  Other  poles,  farther  from  the  origin,  result  from  the 
higher  order  modes  but  the  inverse  transform  of  the  three  dipole  mode 
poles,  i.e.,  the  frequency  dependent  scattering  predicted,  has  been 
shown^^  to  yield  a  fairly  accurate  ramp  response  waveform  for  the  con¬ 
ducting  sphere. 
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Change  in  polarizabiliw  due  lo 
ground  when  dipole  cnodes  do 
not  interact 


If  the  scatterer  is  located  above  or  within  a  half- space  then  an 
interaction  between  the  target  and  the  interface  occurs  which  alters  the 
target  response.  The  simplest  case  to  consider  is  that  of  a  thin  pro¬ 
late  spheroid  which  supports  only  an  electric  dipole  mode  along  the 
major  axis  .  The  poles  of  the  thin  prolate  spheroid  as  a  function  of 
eccentricity  have  been  Cclculated  by  Page  and  Adams,  and  the  electric 
polarisability  has  the  following  form 


^%(8)  = 


{s  +y)^  + 


where  the  complex  poles  are  located  at  s  =  -y  +  jP  (y  and  p  real).  If  an 
electric  field,  Eo(s),  is  applied  parallel  to  the  axis  of  the  spheroid  then 
the  dipole  moment,  Pe(s),  consists  of  a  free-spac#»  term  plus  an  infinite 
sume  of  interaction  terms 


Pe{s)  =  Eo(s)  ae(s)  ^  [ao(s)s  G(s)] "  , 


where  G(s)  is  the  field  scattered  back  from  the  interface  to  a  dipole 
source,  which  in  this  case  is  the  excited  spheroid,  for  a  unit  current 
moment.  The  free-space  term  is  represented  by  n  =  0,  and  remaining 
terms  are  n-th  order  interactions .  Ther  term,  s,  in  Eq.  (48)  is  re¬ 
quired  because  G(s)  is  the  field  scattered  back  for  a  unit  current  moment 
whereas  s  G(s)  is  the  field  scattered  back  for  a  unit  dipole  moment . 
Equation  (48)  can  be  summed  in  closed  form  as 


Pe(®)  =  Eo(s) 


%(s) 

1  -  ae{s)s  G(s) 


Consequently,  the  target  in  the  presence  of  the  interface  can  be  considered 
as  having  a  new  polarizability,  cIq‘{s),  given  by 


“e'(8)  =  T 


a.'s) 


i  =  “c(s)s  G(s) 
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Thus  the  new  polarizability  depends  on  the  target  location  and  orientation 
and  the  half- space  composition,  and  these  effects  are  included  in  C{s]» 

As  the  height  of  the  target  above  the  interface  approaches  infinity,  G(s) 
goes  to  zero  and  ^^'(s)  goes  to  c£g(s). 


3.  Interaction  of  coincident  dipole  modes  - 
in  the  presence  of  the  ground 

b 


The  new  polarizability  given  in  Eq#,(50)  is  valid  for  any  dipole 
mode  which  does  not  interact  with  any  other  dipole  mode*  For  example, 
a  vertical  electric  dipole  produces  only  a  vertical  electric  field  and  no 
magnetic  field  on  the  vertical  axis  and,  as  a  result,  will  not  interact 
with  any  other  dipole  mode*  By  the  same  reasoning,  a  vertical  magnetic 
dipole  mode  will  not  interact  with  other  dipole  modes*  However,  hori¬ 
zontal  electric  and  magnetic  dipole  modes  which  are  perpendicular  to 
each  other  will  interact  because  each  mode  produces  fields  which  excite 
the  other.  The  perpendicular  electric  dipole  moments,  p.,  and  magnetic 
dipole  moments,  p^,  can  be  written  in  the  following  matrix  form 


(51) 
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where  is  the  return  electric  field  produced  by.  a  unit  electric  current 
moment,  is  the  return  electric  field  produced  by  a  unit  magnetic 

current  moment,  Gj^q  is  the  return  magnetic  field  produced  by  a  unit 
electric  current  moment,  and  Ghhi  is  the  return  magnetic  field  produced 
by  a  unit  magnetic  current  moment*  The  matrix  summation  in  Eq*  (51) 
can  be  written  in  closed  form  as 
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It  can  be  seen  from  Eq.  (52)  that  each  dipole  mode  is  influenced  by  both 
the  incident  electric  and  magnetic  fields  in  contrast  with  the  uncoupled 
cases  covered  by  Eqs*  (49)  and  (50), 


oi  New  Dipole  Modes 
in  Target  Detection 

1  •  Use  of  difference  equation 

If  a  target’s  scattering  response  can  be  written  as  a  rational 
function  in  s,  then  the  inverse  transform,  R(t),  must  satisfy  the  follow¬ 
ing  difference  equation  as  shov/n  by  Corrington^^ 

n 

(53)  R(t)=  J  F„,k  R(t-kAt)  , 

k=l 

v/here  the  order,  n,  of  the  difference  equation  is  equal  to  the  number  of 
poles  in  the  Laplace  transform  of  R(t)  and  the  difference  coefficient, 

^n,  k»  simple  functions  of  the  poles  and  the  time  increment,  At. 

The  difference  equation  holds  for  all  values  of  t  >  nAt,  and  the  time 
derivatives  of  R{t)  also  Satisfy  the  same  difference  equation.  Conse¬ 
quently,  the  ramp,  step,  and  impulse  response  waveforms  of  a  target 
each  satisfy  the  same  difference  equation. 

The  exact  solution  for  the  frequency  response  of  the  conducting 
sphere  yields  an  infinite  number  of  poles,  but  it  has  been  shown  that 
the  third  order  difference  equation  resulting  from  the  three  dipole¬ 
mode  poles  given  earlier  quite  accurately  predicts  the  free  space  ramp 
response  of  the  sphere.  If  the  dipole -rneda  poles  of  a  general  scatterer 
in  free  space  are  known:  then  it  is  possibie  to  use  Eqs.  (50)  and  (52) 
to  find  the  new  poles  of  the  scatterer  in  the  preserice  of  ground. 

The  case  of  a  thin  prolate  spheroid  over  a  perfectly  conducting 
ground  is  the  sinipiest  to  ed^^  as  noted  earlier,  because  only  one 

electric  dipble jrnode  is  involyed.  The  polarizability  in  the  direction  of 
the  maior  axis  is  given  in  Eq.  (ir?),"  and  the  new  polarizability  can  be 
found  froni  Eq.  (Sh^^  For  Vertical  orientation  of  the  spheroid  at  a 
height,  h,  abovq^^^t^  as  shown  in  Fig^  11,  G(s)  can  be  shown  to 

have  the  following  form  using  image  theory 

s  +  ” 

(54)  G(s)  =  .- - 

SttCq  c  h^  s 

G(s)  -  =* 

Bit  Cq  c^hs  1.  ; 
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Fig.  H.  Vertical  prolate  spheroid  above 
conducting  ground. 


where  a  power  series  expansion  of  has  been  used  to  obtain  the 

approximate  form  of  G(s).  Using  Eqs .  (47),  (50),  and  (54),  the  new 
polarizability,  cte’(s),  is  written  as 

(55)  a^'(s)  - 


v/here  the  new  poles  of  the  spheroid- ground  configuration  are  the  zeros 
of  the  denominator  in  Eq.  (55),  For  a  spheroid  of  semimajor  axis,  a, 
and  semiminor  axis,  0.013a,  the  free  space  poles  are  located  at 
s  ~  (-  0.1495  +  j  1 , 536)  c/ a.  This  same  spheroid  located  vertically  on 
the  ground  (h=  a),  has  poles  at  s  =  (-0.136ij  1.45)c/a.  The  change  in 
the  poles  due  to  the  presence  of  the  ground  is  not  as  great  as  with  wider 
objects  because  the  very  narrow  spheroid  is  not  a  strong  scatterer  and 
the  interaction  with  the  ground  is  not  strong. 

In  order  to  find  the  new  poles  of  a  conducting  sphere  over  a  ground, 
it  is  necessary  to  consider  both  coupled  and  uncoupled  modes.  The 
vertical  electric  and  magnetic  dipole  modes  are  uncoupled,  and  the  new 
poles  can  be  found  using  the  method  of  Eq.  (50).  However,  the  orthogonal 
horizontal  electric  and  magnetic  dipole  modes  are  coupled  and  the  new 
poles  of  the  configuration  must  be  found  from  Eq.  (52).  Since  the  poles 
result  from. the  inverse  matrix,  they  are  obtained  from  the  zeros  of  the 
determinant: 


SttCq  ch^ 


)  s"  +  2ys  +  +  p*  -  ] 
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(56) 


1  -  CgS  -  0^8  Gg^ 

*“es  ^He  *  -  Gj^  k 

The  solutions  of  Eq.  (56)  for  a  sphere  on  a  conducting  ground  have  been 
obtained  by  again  taking  the  low  frequency  approximation  for  the  G*s  and 
solving  the  resultant  third  degree  equation  in  s.  The  resultant  new  poles 
for  the  coupled  horizontal  modes  for  a  sphere  are  s  =  -  0.898  c/a,  and 
(-  0.801  +  j  0.905)c/a#  The  new  poles  for  the  vertical  electric  mode  of 
the  splp^re  on  a  conductor  are  s  =  (-  0.334±  j  0.626)c/ a,  and  the  new  pole 
for  the  vertical  magnetic  mode  is  s  =  -  0.875  c/a.  Consequently,  the 
sphere-ground  configuration  contains  a  total  of  six  poles  compared  to 
the  three  free- space  poles  given  earlier,  some  of  which  may  not  be 
excited  for  particular  excitation  fields.  For  exainple,  vertical  polari¬ 
zation  of  the  incident  field  will  not  excite  the  vertical  magnetic  mode, 
and  horizontal  polarization  will  not  excite  the  vertical  electric  mode. 

To  check  this  case,  experimental  backscatter  measurements  at  four 
odd  harmonic  frequencies^^  were  used  to  synthesize  the  ramp  response 
of  a  sphere  lying  on  a  conducting  ground.  The  waveform  did  not  fit  the 
free  space  difference  equation  well,  but  it  did  fit  the  free  space  differ¬ 
ence  equation  resulting  from  the  new  pOles  quite  well,  particularly  for 
vertical  polarization. 

If  the  difference  equation  for  a  particular  object  is  known,  it  can 
provide  the  following  detection  scheme.  The  received  waveform  is 
sampled  in  time  and  the  difference  equation  of  the  desired  object  is 
used  to  predict  a  new  waveform.  If  this  waveform  agrees  closely  with 
the  received  waveform,  then  the  presence  of  the  desired  object  is  in¬ 
dicated 

Note  particularly  that  an  analytical  sohition  for  the  poles  of  the 
target  on  the  ground  is  not  a  prerequisite  for  utilization  of  this  pre¬ 
dictor  type  processing.  The  method  requires  only  that  the  difference 
equation  for  the  ramp  response  waveform  of  a  target  on  the  ground  be 
known.  The  constants  in  the  difference  equation  may  be  deter¬ 

mined  from  an  analysis  of  experimental  multiple  frequency  scattering 
data  on  a  particular  target.  A  rational  function  approximation  to  the 
complex  scattering  amplitude  as  a  function  of  frequency  is  formed, 
wherein  the  order  n  of  the  difference  equation  is  equal  to  the  number 
of  poles  in  the  rational  approximate  and  the  coefficients  F^^^  k  are  simple 
functions  of  the  pole  locations  and  the  time  increments,  At.  For  spherical 
scatterers  and  a  few  other  shapes,  the  pole  locations  can  be  determined 
analytically,  and  a  rational  approximant  of  any  order  obtained. 
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The  basic  predictor  or  difference  equation  technique  was  tested 
by  trying  to  detect  a  hemispheric 2I  boss  on  a  conducting  ground  in  the 
presence  of  artificial  trees.  Backscatter  measurements  at  four  odd 
harmonics were  used  to  synthesize  a  ramp  response  waveform,  and 
the  difference  equation  for  the  hemispherical  boss,  which  is  the  same 
as  the  difference  equation  for  the  free  space  sphere,  was  used  to  predict 
a  new  waveform  •  The  correlation  of  the  predicted  and  received  wave¬ 
forms  for  vertical  polarization  was  higher  for  the  trees  plus  the  boss 
than  for  the  trec»  alone  even  for  a  fairly  xow  signal  to  clutter  ratio. 

For  example,  using  a  time  increment  of  0.8  c/a,  the  addition  of  the  boss 
to  the  trees  produced  an  increase  in  correlation  from  0.382  to  0.734 
when  the  signal  to  clutter  power  ratio  was  1.72.  The  increase  in  cor¬ 
relation  was  greater  for  higher  signal  to  clutter  ratios  and  smaller  for 
lower  signal  to  clutter  ratios.  Consequently,  a  proper  threshold  setting 
for  the  correlation  would  permit  detection  of  the  boss  in  the  trees.  One 
advantage  of  this  method  is  that  even  though  the  backsc altered  field  is  a 
function  of  aspect  angle,  the  difference  equation  is  not.  Consequently, 
an  integration  on  the  aspect  angle  could  be  accomplished  without  changing 
the  difference  equation  receiver. 

To  test  the  discrimination  properties  of  the  difference  equation, 
the  theoretica.1  free- space  ranip  responses  of  Spheres  either  larger  or 
smaller  then  the  desired  size  sphere  were  fed  into  the  receiver.  Using 
a  At  of  1.6  c/a,  the  correlation  dropped  from  1.0  to  0.5  when  the  sphere 
radius  was  increased  55%  or  decreased  20%.  The  correlation  continues 
to  drop  as  the  sphere  size  is  changed  more,  but  the  drop  is  more  rapid 
with  decreasing  size.  This  is  expected  because  it  has  been  generally 
observed  that  rapidly  varying  waveforms  are  easier  to  discriminate 
against  than  slow^ly  varying  waveforms. 

2.  Backscattered  field 

The  presence  of  the  ground  alters  the  low  frequency  scattering 
properties  of  an  object  by  modifying  the  polarizability  of  the  dipole  modes 
and  by  changing  the  incident  field.  The  first  effect  alone  detei*mines  the 
shift  in  the  poles  and  is  sufficient  for  difference  equation  calculations 
which  involve  only  the  poles  and  not  the  excitation.  However,  if  the  back- 
scattered  field  as  a  function  of  frequency  and  aspect  are  required  as  in 
a  matched-filter  receiver,  then  both  effects  are  important. 

To  illustrate  a  method  of  combining  the  two  effects,  consider  the 
case  of  a  thin  prolate  spheroid  vertically  oriented  at  a  height  h  above 
ground  (Fig.  11).  For  a  vertically  polarized  plane  wave,  the  vertical 
component  of  the  incident  field  at  the  center  of  the  splieroid  is  given  by 


2h 

-s  —  cos  6 

£2(5)  =  £q(8)  1  +  Ry  ® 


sinO 


where  0  is  the  aspect  angle  measured  from  the  vertical  and  Ry  is  the 
reflection  coefficient  of  the  ground  for  vertical  polarization.  Using 
Ekjs.  (50)  and  (57),  the  dipole  moment  of  the  spheroid  is 


p(s)  = 


l-a(s)s  G(s) 


The  far  zone  backscattered  field  in  terms  of  the  dipole  moment  is 
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Using  Eqs.  (58)  and  (59)  the  complex  pha&or  response,  G®(s),  of  the 
spheroid  over  the  ground  is 


sin^O 


2h 

*s  — cos  0 
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°  ■  4Tr€oC='  l.a(s)sG(s) 


The  two  effects  of  the  ground  on  G®(s)  can  be  seen  directly  from  Eq.  (60), 

The  change  in  polarizability  is  represented  by  the  factor,  [1  -  a(s)s  G(s)] 
and  the  altered  excitation  is  represented  by  the  factor,  [  1  -S- Ry  e"®  ^h/ c  cosojz^ 
The  problem  is  more  complicated  for  objects  v^hich  may  require  additional 
dipole  modes,  but  the  procedure  can  be  extended  using  the  modified 
polarizabilities  in  Eqs.  (50)  and  (52). 


VII.  CONCLUSIONS 


Exact,  closed  form  solutions  for  the  transient  fields  of  vertical 
and  horizontal  electric  dipoles  in  the  presence  of  a  half- space  have  been 
obtained  for  the  special  case  of  a  nondispersive  half- space  and  particular 
source-observer  orientations.  The  same  technique  used  to  obtain  these 
results  can  be  used  to  obtain  exact  closed  form  solutions  for  vertical  and 
horizc  .tal  magnetic  dipoles  under  the  same  conditions,  but  these  solutions 
have  not  been  carried  out  in  this  report.  For  a  nondispersive  half- space 
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but  an  arbitrary  observer  location,  the  formulation  leads  to  a  finite 
integral  which  must  be  evaluated  numerically. 

If  the  half-space  is  dispersive  or  btratifieds  exact  closed  form 
solutions  are  not  possible  for  the  transient  fields  and  approximate 
techniques  must  be  utilized*  For  a  dispersive  half-space  and  observa¬ 
tion  points  above  on  the  half- space,  the  form  of  the  transient  fields 
is  such  that  reasonable  estimates  are  possible  from  a  combination  of 
short  and  long  time  approximations*  Some  of  the  necessary  approxi¬ 
mations  have  been  obtained  using  both  quasi- static  and  asymptotic 
approximations  and  a  Fourier  synthesis  technique. 

A  general  model  for  the  backscattering  including  interactions 
from  electrically  small  (Rayleigh  and  low  resonance  range)  conducting 
targets  located  on  the  ground  has  been  formulated.  The  model  illustrates 
that  the  lov/  frequency  effects  of  the  ground  on  the  scattering  properties 
of  the  target  are  a  modification  of  the  polarizability  of  the  dipole  modes 
of  the  target  and  an  alteration  of  the  incident  field.  These  effects  have 
been  illustrated  for  the  case  of  a  thin,  prolate  spheroid  target. 

A  predictor  type  difference  equation  has  been  suggested  for  target 
detection  and  discrimination.  The  method  requires  that  the  difference 
equation  for  the  ramp  response  of  a  target,  either  in  free  space  or  on 
the  ground,  be  known  and  the  equation  can  be  determined  from  multiple 
frequency  scattering  data.  For  a  few  target  shapes,  the  difference 
equation  can  be  determined  a  lalytically.  Application  of  the  difference 
equation  for  detection  has  been  illustrated  using  simulated  measured 
data  of  conducting  hemispherical  boss  targets  on  the  ground  under-  trees. 
Discrimination  a.ppliccition  has  been  illustrated  for  size  separation  of 
free  space  spherical  scatterers. 

Research  on  the  problems  discussed  in  this  report  is  continuing  on 
Contract  FI  9628- 67 -  C-0239.  A  detailed  report  including  extended  re¬ 
sults  for  both  the  transient  fields  of  dipoles  in  the  presence  of  a  half¬ 
space  and  applications  of  the  difference  equation  to  detection  and  dis¬ 
crimination  is  anticipated  in  the  near  future. 
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